Heating and atom loss during upward ramps of Feshbach resonance levels in Bose-Einstein 

condensates 
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The production of pairs of fast atoms leads to a pronounced loss of atoms during upward ramps of Feshbach 
resonance levels in dilute Bose-Einstein condensates. We provide comparative studies on the formation of 
these bursts of atoms containing the physical predictions of several theoretical approaches at different levels of 
approximation. We show that despite their very different description of the microscopic binary physics during 
the passage of a Feshbach resonance, all approaches lead to virtually the same prediction on the total loss of 
condensate atoms, provided that the ramp of the magnetic field strength is purely linear. We give the reasons 
for this remarkable insensitivity of the remnant condensate fraction to the microscopic physical processes and 
compare the theoretical predictions with recent Feshbach resonance crossing experiments on 23 Na and 85 Rb. 

PACS numbers: 03.75.Kk, 34.50.-s, 05.30.-d 



I. INTRODUCTION 

Since their first experimental demonstration |.l], magnetic 
field tunable interactions have proven their potential and use- 
fulness in the physics of quantum degenerate gases. This 
technique takes advantage of the Zeeman effect in the elec- 
tronic energy levels of the atoms to manipulate their binary 
collision properties in the vicinity of a Feshbach resonance. 
For several atomic species it is now possible to change the 
scattering length by orders of magnitude, and even to reverse 
its sign, making the atoms either repel or attract one another. 
One of the impressive applications of this technique was the 
successful creation of stable 85 Rb Bose-Einstein condensates 
with thousands of atoms |2], that would otherwise collapse if 
consisting of more than about 80 atoms. These experiments 
have now been refined to associate pairs of atoms to ultracold 
molecules a nd p robe their coheren ce prop erties 010,1110,0 
II S El [H El El El El El El E USUI. The dissocia- 
tion energies of the highly excited dimer molecules have been 
measured accurately down to several kHz HHH). 

Early attempts to associate Feshbach molecules in Bose- 
Einstein condensates 1 2 1 ] relied upon linear ramps of the mag- 
netic field strength across the resonance position, in which the 
resonance level (see, e.g., I22l 12311 ') was s wep t upward in en- 
ergy 12411 . Subsequent theoretical studies 1 22, 25] have shown 
that these particular ramps lead to a heating of the cloud rather 
than to molecular production. The condensate atoms were lost 
into correlated pairs of atoms with a comparatively high veloc- 
ity of their relative motion. These atom pairs were also pro- 
duced in the interferometric experiments of Ref. (3[, involving 
a special sequence of magnetic field pulses, and directly de- 
tected as a dilute burst of atoms with a spread far beyond the 
size of the remnant Bose-Einstein condensate. The theoretical 
studies in Refs. 126. 1271 12811 have confirmed this interpreta- 
tion and provided a quantitative understanding of the number 
of burst atoms observed in Ref. H] as well as their measured 
energy spectrum 1 28[ . 

Several approaches have been proposed to describe the con- 



version of condensate atoms into bound or continuum molec- 
ular states (see, e.g., m m m m m m m m m m 

I33II34I.I35I0 . Early approaches are based either on mean field 
theory, with an exchange of particles between the atomic con- 
densate and the Feshbach resonance level (e.g. 1 251 l29l.l30ln . 
or on a two-body configuration interaction picture, extended 
in such a way that it includes the Bose stimulation of the tran- 
sition between free atom pairs and the resonance level 1 22] . 
While these approaches effectively include only two energy 
levels of the atom pairs, i.e. the zero energy scattering state 
and the Feshbach resonance level l22T.l23ll . the interferometric 
experiments in Ref. 1 3J have clearly revealed that a complete 
description of the molecular association requires many-body 
approaches, which treat the two-body collision dynamics in- 
cluding all states in a non-perturbative manner j^lSIH]- 

In this pap er we use the beyond mean field approach of 
Refs. 128113211 to describe the atom loss in upward ramps of a 
Feshbach resonance level in a Bose-Einstein condensate and 
compare the results with the simpler two level mean field ap- 
proach of Refs. I25II29II3O I. Our studies involve the Feshbach 
resonance crossing experiments of Refs. |21] and |2]. In the 
experiments 1 2 1 ] a Bose-Einstein condensate of 23 Na in the 
(F = l,m F = +1) hyperfine state was exposed to a mag- 
netic field swept across the Feshbach resonances at 853 and 
907 G 1 1], respectively. In Ref. |2] the magnetic field strength 
was swept across the 155 G resonance of a gas of 85 Rb atoms, 
which was condensed in the (F = 2, nif = -2) hyperfine state. 
In all cases dramatic losses of atoms were observed, which in- 
creased with decreasing ramp speed. 

The widths of the 23 Na Feshbach resonances have been pre- 
dicted to be 1 G (at 907 G) and 0.01 G (at 853 G), respectively 
E|. The 85 Rb 155 G Feshbach resonance has been charac- 
terised very accurately in several experimental and theoretical 
studies I3ll4l l3a.l37ll . and its width is predicted to be 1 1 G. We 
shall show on the basis of binary collision physics that the mi- 
croscopic dynamics caused by the ramps of the magnetic field 
strength is decisively different between the 23 Na experiments 
l2lll and the experiments |2] with the broad 85 Rb Feshbach 
resonance, although the direction of these resonance crossings 



2 



excludes the adiabatic association of atoms to molecules (see, 
e.g., I22l 12311 ) in both cases. We shall show that a two level 
picture is particularly inappropriate to describe the Feshbach 
resonance crossing in 85 Rb |2]. Our studies reveal, however, 
that both the first order microscopic quantum dynamics ap- 
proach of Re fs. |28[l32ll and the two level mean field theory of 
Refs. 0H 30] predict largely the same atom loss for the ex- 
periments in Refs. 12111 and |2]. We show that this remarkable 
insensitivity of the remnant condensate fraction with respect 
to the description of the microscopic binary collision physics 
is characteristic for linear ramps of the magnetic field strength 
and provide the essential physical parameters that determine 
the losses. We show furthermore that our findings are quite 
general as they also apply to the dissociation spectra of Fes- 
hbach molecules in upward ramps of the Feshbach resonance 
level, which are widely used for the purpose of molecular de- 



Here H™(B) = -H 2 V 2 /2m + V* ap + E™(B) is the Hamilto- 
nian of a single atom in the hyperfine state jj containing the 
kinetic energy, the trap potential and the hyperfine energy that 
depends on the magnetic field strength B due to the Zeeman 
effect, and m is the atomic mass. V^vj.j^K 1 ") is the potential 
associated with the asymptotic incoming and outgoing binary 
scattering channels \kA) and \jiv], which depends on the rel- 
ative coordinates r = x - y. The field operators obey the 
Bose commutation relations [(/^(x), (jj\,{y)] = 5^ v 5(x - y) and 
[^(x),^ v (y)] = 0. 

All physical quantities of a gas are determined by cor- 
relation functions, i.e., expectation values (denoted by 



Here H™ ][kA) {B) = [H l *{B) + H v B {B)]6 KIA 6 Av + is 
the Hamiltonian matrix of two atoms. All the other correla- 
tion functions can be expressed in terms of the mean field and 
the pair function |28]. The coupled equations (0 and Q are 
the most general form of the first order microscopic quantum 
dynamics approach of Refs. 1281 13211 . but they also include 
the two level mean field approach of Refs. 1251 l29i 13011 in 
the Markov approximation |23]. Higher order approximations 
lEoll in the general truncation scheme l32ll include few-body 



tection (see, e.g., BSE l38ln . We note, however, that the 
interferometric experiments of Refs. |3j,|4J] provide an example 
of a sequence of magnetic field pulses, which as a whole varies 
non-linearly in time. These experiments indeed reveal details 
of the microscopic binary collision physics beyond universal 
considerations MEslIm 



II. MICROSCOPIC QUANTUM DYNAMICS APPROACH 

In general, the interatomic interaction depends on the rel- 
ative coordinates of the atoms and couples all their hyperfine 
energy levels, labelled in the following by Greek indices. A 
trapped gas of Bose atoms is then described by the many -body 
Hamiltonian: 



(1) 

I 

(••■)) of normal ordered products of field operators in 
the state at time t. Quantities of particular interest are 
the mean field ^(x, f) = ((/^(x)),, the pair function 
% v (x, y, t) = (if/yiy^ix)),-^^, f)^(x, t) and the one-body 
density matrix of the non-condensed fraction r jUV (x, y, f) = 

<^(y¥/Xx)>, - T;(y, tyV^X, t). We have shown in Refs. 
13211 how the exact infinite hierarchy of dynamic equations for 
correlation functions can be truncated to any desired degree 
of accuracy. In the first order approach 1 28] the mean field 
and the pair function are determined by the coupled dynamic 
equations l; 23ll : 



(2) 
(3) 



inelastic loss phenomena, like three-body recombination |41 ]. 
We shall focus in this work on the dominant two-body scatter- 
ing in the presence of the surrounding Bose gas. 



H = 2 f dx 0J«#™(*WmGO + f dx f dy <A;(x¥j(y)Vw,M)(x - y)«x)Wy). 



I 

ih^V^x, t) =//; B (B)^(x, t) + J] [dy %(y, f)V w , M (x - y) [O^x, y, t) + ^(x, t)V A (y, t)] , 

vkA J 

ih-% v (x, y, t) = 2 KW B )°^( X - y. f ) + *W*>(* - y)^( x > f )^(y, ■ 

kA 

I 
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III. RESONANCE ENHANCED BINARY COLLISIONS 

A. General properties of Feshbach resonances in ultracold 
collisions 

We denote the open s wave binary scattering channel of 
two asymptotically free 23 Na condensate atoms in the {F = 
\,Mf = 1) state as the {1,1} open channel (or simply as the 
open channel) with an associated background scattering po- 
tential V(i,i},{i,i}(r) = Vhg{r), where r is the internuclear dis- 
tance. Analogously, for 85 Rb the background scattering po- 
tential is denoted by V(-2,-2),(-2 ; -2)( r ) = Vb g (r). Throughout 
this paper we shall choose the zero of energy as the dissocia- 
tion threshold of VbgM, i-e. VWW ~ 0. Figuref^shows, for 

the case of the 85 Rb Feshbach resonance, that in the vicin- 
ity of the range of magnetic field strengths (from 162.2 to 
132 G) in Ref. |2] the energy E KS (B) of a closed channel 
vibrational state (a Feshbach resonance level) res ('") crosses 
zero (at B res = 164.1 G). The resulting strong coupling be- 
tween the open and the closed channel is accompanied by 
the emergence of a highly excited multichannel vibrational 
molecular bound state 0b, whose binding energy E\,{B) van- 
ishes at the measurable position of the Feshbach resonance 
of Bo = 154.9 G. When the magnetic field strength is swept 
across the resonance from above, i.e. the energy of the Fesh- 
bach resonance level <p res (r) is tuned upward, 0b thus ceases to 
exist (see Fig.Q. The s wave scattering length of two asymp- 
totically free condensate atoms, therefore, has a singularity at 
Bo. The emergence of a multichannel diatomic bound state, 
which causes the singularity of the scattering length at Bo, is 
the decisive feature of all Feshbach resonances in ultracold 
collision physics. 



B. Two channel approach 

The coupling between the {-2, -2} open channel of 85 Rb 
and other open channels, leading, e.g., to inelastic two-body 
losses, is weak |41]. In the case of 23 Na the condensate atoms 
are in the electronic ground state and inelastic two-body losses 
are ruled out from the outset. It thus turns out that for the 
Feshbach resonances of both species the Hamiltonian of the 
relative motion of two atoms in free space can be reduced to 
just two channels J22ll43ll : 



-v 2 [|bg)<bgt 

m 



H 2B {B) 
V 2B (B)= (|bg>,|cl» 



cl>(cl|] + 

V be (r) W(r) 
W(r) V cl (B,r) 



V 2 b(B), (4) 

<bgl ) (5) 
<cl| r { > 



Here |bg) is the product of hyperfine states associated with the 
open channel, |cl) is the superposition of products of atomic 
hyperfine states associated with the strongly coupled closed 
channel [22], and W(r) couples the two channels. V C \(B, r) is 
the closed channel potential that supports the resonance state 
4> KS (r), i.e. 



[-H 2 V 2 /m + V c i(fl,r)]0 res (r) = E KS (B)^ res 



sir). (6) 
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FIG. 1: Magnetic field dependence of the binding energy Ey,(fi) 
of the highest excited multichannel vibrational bound state <p b of 
85 Rb2 (solid curve) and the energy E ICS (B) of the closed channel res- 
onance state (f) ICS (r) (dashed curve). £ rcs (6) depends linearly on B 
with fr l dE as /dB = -3.46 MHz/G. The inset shows E b (B) (solid 
curve) as compared to the universal asymptotic formula -ft 2 /ma 2 for 
the near resonant binding energy (dotted curve) and its first improve- 
ment l42ll that accounts for the length scale set by the van der Waals 
interaction in a single channel picture (dotted dashed curve). 



As only the resonance state couples strongly to the open chan- 
nel, -n 2 V 2 /m + V c i(B,r) can be replaced by the one dimen- 
sional Hamiltonian |0 res )£Ves(BX0resl to an excellent approx- 
imation l23ll . Within the limited range of magnetic field 
strengths under consideration the energy associated with the 
resonance state is linear in the magnetic field strength, i.e. 



Eres(B) 



dF 



dB 



-(■Bres) 



(B - fires) ■ 



(7) 



Here dE les /dB denotes the difference between the magnetic 
moment of the resonance state and an atom pair in the open 
channel, and is usually referred to as the slope of the reso- 
nance. 



C. Effective low energy potentials 

The typical de Broglie wave lengths of the condensate 
atoms are much larger than the length scales set by the spa- 
tial extent of the binary interactions. The details of the po- 
tentials are, therefore, not resolved in the binary collisions. 
Consequently, we can use the approach of Ref. 12311 whose pa- 
rameters for the 23 Na and 85 Rb Feshbach resonances we shall 
briefly summarise. For the potential energy operator Vb g we 
use the separable form 



v bg = hfbg>& g <Ybgl> 



(8) 
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with the convenient Gaussian ansatz 



( V2^o- bg ) 3 ' 



(9) 



as introduced in Ref. 12 811 to describe the low energy back- 
ground scattering of 23 Na and 85 Rb SHUH^l- The 
background scattering potential Vb g then recovers the back- 
ground scattering length of the atoms in the vicinity of the 
respective resonance, as well as the length scale 



vdW 



(10) 



set by the van der Waals interaction -C(,/r 6 , which is de- 
termined by the dispersion coefficient Ce- Table H] provides 
the relevant physical parameters associated with the Feshbach 
resonances considered in this paper, and Table shows the 
parameters an d °"bg characterising the low energy back- 
ground scattering potential. To recover the widths AB of the 
resonances and the shifts Bo - Z? res (see Fig.QJ as given in Ta- 
ble|l]it is sufficient to specify W(r)<p res (r). We have used the 
ansatz 



W(r)0 res (r) = ( 



( V2^o-) 3 ' 



(11) 



The parameters £ and <x can be obtained from Table [H] The 
scattering length is then exactly given by the formula 



a(B) - abg 1 - 



AB 
B-B 



(12) 



The binding energies E\,(B) in Fig.[2have been obtained with 
this low energy potential matrix V2b(B) by 



H2B<Pb = Eb<f>b- 



(13) 



Figure |3 shows a comparison of E\,(B), for the 907 G Fes- 
hbach resonance of 23 Na, between the low energy poten- 
tials in this work and the highly accurate two channel model 
of Ref. l22ll . The associated two-body time evolution op- 
erator U2B{t,T) can be obtained by relatively simple means 
ll28ll . even for the time dependent magnetic field strengths in 
Refs. 1I2 H21I1 . U2B(t, t) will serve as an input to the solution of 
the coupled equations and 0. 



D. Differences between the 23 Na and 85 Rb Feshbach 
resonances 

Although the 23 Na and 85 Rb Feshbach resonances consid- 
ered in this paper can all be described by a low energy two 
channel Hamiltonian of the form of Eq. (0}, with a single res- 
onance state, the multichannel bound states which cause the 
singularities of the scattering length differ decisively among 
the resonances within comparable ranges of the magnetic field 
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FIG. 2: Magnetic field dependence of the modulus of the binding 
energy of the highest excited multichannel vibrational bound state of 
23 Na2 in the vicinity of the 907 G Feshbach resonance. The fig- 
ure compares the binding energies supported by the two channel 
low energy potential matrix of Ref. 1 23] and this work (solid curve) 
with the highly accurate two channel model of Ref. 12211 (filled cir- 
cles). The dotted curve indicates the universal asymptotic formula 
E\, = -h 2 /ma 2 and the dashed curve is the Gribakin and Flambaum 
l42ll estimate that accounts for the length scale set by the van der 
Waals interaction in a single channel picture. We note that the bind- 
ing energies are given on a logarithmic scale. 



strength. The differences are illustrated in Fig.|5J which shows 
the admixture |(0 res , cl|0b)| 2 of the closed channel state to 
the two channel vibrational molecular bound state 0b of the 
Hamiltonian @. For 85 Rb the admixture stays small (< 30 
%) over the entire range of experimental |2] magnetic field 
strengths from 154.9 to 162.2 G on the high field side of the 
Feshbach resonance (cf., also, Ref. |43]). Even for the broader 
(AB = 1 G) resonance of 23 Na at 907 G, however, the mul- 
tichannel bound state 0b is transferred immediately into the 
closed channel resonance state <p res (r), with an admixture of 
more than 80 % within the resonance width. The single chan- 
nel nature of the 85 Rb 155 G Feshbach resonance is also visi- 
ble in the binding energies in Fig.^because the near resonant 
single channel Gribakin and Flambaum estimate l28ll42Tl 



£b — —ti 2 lm(a — a) 2 



is applicable, where 



a - 



1 rp/4) 

V2H5/4) 



(14) 



(15) 



is the mean scattering length |42], which accounts for the long 
range behaviour -Ce/r 6 of the background scattering poten- 
tial in terms of the van der Waals length in Eq. dlOl l. For the 
907 G Feshbach resonance of 23 Na, however, Fig.|5]shows that 
Eq. dl4> is even less accurate than the universal near resonant 
estimate j3H3E3] £t> = -h 2 /ma 2 of the binding energy. 

The reason for the single channel nature of the 85 Rb 155 G 
Feshbach resonance is the large shift of Bq - B res = -9.2 G 
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Element a bg [a B ] C 6 [a.u.] [MHz A] fi [G] AS [G] B - B res [G] dE res /dB [MHzh/G] 

23 Na 63.9 [22] 1561 [44] 207.8 [45J 853 [21J 0.01 [22] -5.85 xlO" 3 [46J 5.24 [22] 

62.8 [22] 1561 [44] 207.8 [45J 907 [21] 1.0 [22] -0.55 [46] 5.24 [22] 

85 Rb -450 [3] 4660 [36J 155 T36] 11.0 [36J -9.2 [46] -3.46 [22J 

TABLE I: Parameters characterising the Feshbach resonances (aaohr = 0.052918 nm; Q: 1 a.u.= 0.095734 x 10~ 24 J nm 6 ). 

Element B [G] £ bg [4rffl Bohl /m] cr bg [a B ohr] f [GHz ftaj^J cr [a Bohr ] 

23 Na 853 ^197 27^2 ^96 14.5 

907 -186 26.5 -44.0 17.4 

85 Rb 155 -99.3 71.9 19.0 52.4 

TABLE II: Parameters determining the low energy two channel potential matrix. 



(see Fig.^, which leads to a large detuning of E tss (Bo)/h = 32 
MHz at the resonance position. The naive physical two level 
picture, in which the resonance level <p K s(r) is swept across 
the zero energy level of the background scattering potential 
at the measurable position Bq of the Feshbach resonance, and 
is then transferred to a bound molecular state (f>\, at negative 
energies E les < (see, e.g., Ref. |i4g]), is therefore inappli- 
cable. In fact, the resonance state does not even cross the 
threshold within the entire range of experimental magnetic 
field strengths (from 162.2 to 132 G) in Ref. Q]. 



that Eq. dl6l accounts for the entire two-body collision dy- 
namics including the continuum of modes above the disso- 
ciation threshold of the background scattering potential in a 
non-perturbative manner. 



B. Contributions beyond mean field approaches 

The solution of Eq. d 1 6t > determines the pair function, via 
Eq. Q, through its components in the open and closed chan- 
nels: 



IV. MANY-BODY APPROACHES TO FESHBACH 
RESONANCE CROSSING EXPERIMENTS 

A. Non-Markovian nonlinear Schrodinger equation 

Inserting the two-channel Hamiltonian @ into Eq. im- 
plies that the pair function has a closed channel component 
3>ci(x, y, t), proportional to <f> res (r), and a long range compo- 
nent C>bg(x, y, f) in the open channel. Only the (F = 2, nif 

1 ) hyperfine state of 85 Rb and the (F = I, trip = 1) hyper- 
fine state of 23 Na have a mean field that we will denote simply 
by ^(x, t). The dynamic equation Q can be solved formally 
in terms of the mean field and inserted into Eq. l28 H32ll . 
In the first order microscopic quantum dynamics approach the 
entire dynamics of the gas is then described by a single non- 
linear Schrodinger equation 



|<D(f)> = |bg><Dhg(f) + |cl)(D cl (f). 



(18) 



i%— W(x,f) 
dt ' 



n 2 7 
- r v 2 + yx) 

2m 

-*P*(x,f) I dr? 

J to 



W(x, t) 



(X,T) — h(t,T), (16) 
OT 



with 



h{t,r) = (27r/0mbg|V 2B (B(r)) f/ 2 B(f,r)|O,bg)0(f- t). (17) 

Here Uzs (t,r) is the time evolution operator associated with 
the magnetic field dependent Hamiltonian @, 9(t - r) is the 
step function, which yields unity for t > r and vanishes else- 
where, and |0, bg) is the zero momentum plane wave of the 
relative motion of two atoms in the open channel. We note 



A derivation |23, 28] beyond the scope of this paper shows 
that the pair function determines the components of the den- 
sity matrix of the non-condensed fraction in the open and 
closed channel to be 

r bg (x,x',f) = Jjy0) bg (x,y,f)O; g (x',y,f) (19) 



and 



r c i(x,x',f) = J dy® ci (x,y,tW cl (x',y,t), 
respectively. The density of the gas is thus given by 

«(x, o = r bg (x, x, t) + r d (x, x, t) + ^(x, oi 2 , 



(20) 



(21) 



and the total number of atoms N - J dx n(x, t) is a constant 
of motion I23l 12811 . 

In general, the solutions of the coupled equations (|2ji and 
Q depend not only on the initial condensate wave function 
^(x, fo) and but also on the initial pair function O(fo), which 
has so far been neglected in Eq. d!6l >. These initial conditions 
are determined by the experimental situation. For Feshbach 
resonance crossing experiments with 23 Na it appears reason- 
able to require the condensate density to be stationary in the 
absence of any variation of the magnetic field strength. To find 
the stationary solution of the coupled equations (|2j and (13, we 
assume that the magnetic field strength at time fo is sufficiently 
far away from the Feshbach resonance that the binary dynam- 
ics is determined entirely by the background scattering. The 
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FIG. 3: Magnetic field dependence of the population of the closed 
channel component of the highest excited vibrational bound state 
<p h associated with the 23 Na 907 G and the 85 Rb 155 G Feshbach 
resonance. The populations were determined from the two channel 
Hamiltonian J4j. 



closed channel component of the initial pair function can then 
be neglected, and for its open channel component <t>b g as well 
as for the mean field *F we choose the ansatz: 



¥(R, t) = V"o(RK i/Jo( '~' o)/fi , 

0> be (R, r, t) = O bg (R, r, fo^ - ^ - ' ^. 



(22) 
(23) 



Here R = (x + y)/2 and r = y - x are centre of mass and rela- 
tive coordinates, respectively, no is the equilibrium condensate 
density, and fiQ can be interpreted as the chemical potential. 
The formal solution of Eq. then provides the general form 
of the stationary pair function 



<D bg (R, r, to) = « (R)(2^) 3/2 Re [<r|G bg (2 W) + i0) V bg |0>] ■ 



(24) 



where 



G bg (2// + i0) = [2fxa + i0 + h 2 V 2 /m - V bg ]~' (25) 



denotes the energy dependent Green's function associated 
with the background scattering potential V bg . Inserting the 
stationary ansatz of Eq. d22b for the mean field into Eq. @ 
then yields the stationary nonlinear Schrodinger equation: 



A'o 



V«o(x) = 



h 2 , 
~2^ Vtrap(x) 



V«o(x) 



[«o(x)] 3/2 (2^) 3 Re [<0|T bg (2 / z + /0)|0>] . 



(26) 



Here r bg (2/in + i0) is the energy dependent T matrix associ- 
ated with the background scattering potential V bg (see, e.g., 
Ref. 1331 ). Evaluating the T matrix at the energy 2/in — > 0, 
Eq. ( I26> recovers the usual Gross-Pitaevskii equation in the 
dilute gas limit (w p ka b <k 1, where w p k is the peak density), 
which then determines the chemical potential [Xq as well as 
the stationary condensate density «o(x). In the experiments 
with 23 Na 1 21] the peak condensate density is on the order of 
n p k ~ 1 .5 x 10 15 cm -3 and including the initial pair function in- 
fluences the dynamics on the time scales of the passage across 
the Feshbach resonances. We have thus included the initial 
pair function in Eq. d!6i as explained in Ref. 1 32] . In the ex- 
periments with 85 Rb the peak density n p k ~ 2.2 x 10 12 cm~ 3 
turns out to be sufficiently low for initial pair correlations to 
be negligible. 



C. Two level mean field approach 

The two level mean field approach of Refs. ESUHH can 
be derived from the first order microscopic quantum dynamics 
approach in Eqs. (|2) and (0 with the ansatz 



cD d (R,r,f)= V2> res (ryF res (R,f) 



(27) 



for the closed channel component of the pair function in centre 
of mass and relative coordinates. HerejA re s(r) is the closed 
channel resonance level. A derivation 1 23] beyond the scope 
of this paper then shows that Eqs. (0 and yield the coupled 
equations 



m-y(x,t) = 

at 



V 2 + — ojL 

2m 2 ho 



W(x, t) 



+ g bg |¥(x,f)| 2> P(x,f) 

+ g KS V*(x,t)^ KS (xj) 

for the condensate mean field, and 



(28) 



iti^ KS (R,t) 
ot 



fl 2 r, 2m , n 



dE K 



dB 



(Bres) 



¥ res (R0 

[B(t)-B ]V Ies (R,t) 



+ - F g«V 2 (R,t) 
V2 



(29) 



for the mean field associated with the resonance level in the 
Markov approximation. Here we have assumed that the atom 
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trap is spherical and that the closed channel resonance state 
experiences twice the trap potential of the atoms. The cou- 
pling constants gbg and g res are given by gbg = Anf^a^/m 
and g res = (27r^) 3 ^ 2 (<^> res |W|0[ ) +) ), where «b g is the background 
scattering length, and (p[^\r) is the zero energy wave func- 
tion associated with the background scattering 1 23 ] . The 
physical picture underlying the two level mean field approach 
(see, e.g., Ref. [46]) can be described in terms of a linear 
curve crossing between the energy E res (B) associated with 
the closed channel resonance level <f> ies (r) and the zero energy 
level of two condensate atoms (see Fig.[0. This picture, how- 
ever, neglects all the other modes of the two-body background 
scattering continuum and, consequently, does not recover the 
nature of the highest excited vibrational diatomic bound state 
lE^ll . Even more importantly, in the context of this paper, the 
absence of continuum modes also implies that the two level 
mean field approach is unsuited to describe the energy dis- 
tribution of the atoms lost from the condensate in an upward 
ramp of the Feshbach resonance level. 

V. COMPARISON BETWEEN DIFFERENT MANY-BODY 
APPROACHES 

A. Loss of condensate atoms in upward ramps of the 
resonance level 

Since, in accordance with Eqs. ( 1191 and J20i . the non- 
condensate density stems directly from the pair function, the 
build-up of pair correlations is the dominant source of atom 
loss from the condensate. In general, a time dependent mag- 
netic field can transfer pairs of trapped condensate atoms ei- 
ther into highly excited diatomic molecules or into the com- 
plementary part of the two-body energy spectrum that con- 
tains a quasi continuum of excited levels 12 811 . In the exper- 
iments in Refs. |2j 12111 molecular formation can be ruled out 
because the linear ramps start on the side of the Feshbach res- 
onance where the scattering length is positive and cross the 
resonance in such a way that the highest excited multichannel 
vibrational bound state <pb ceases to exist during the passage 
(see Fig. [0 1 4-8 ] . The remaining source of atom loss from 
the condensate is the formation of correlated pairs of atoms in 
initially unoccupied excited energy levels. As the time de- 
pendent homogeneous magnetic field provides only energy 
but no momentum the centre of mass motion of these pairs is 
limited by the initial momentum spread of the Bose-Einstein 
condensed gas. Their comparatively fast relative motion, how- 
ever, has been observed as a very dilute burst of atoms, e.g., in 
Ref. and the typical scale of experimental single particle 
kinetic energies of /ikin/^B ~ 150 nK has been confirmed with 
the first order microscopic quantum dynamics approach of this 
work 12811 . We note that the two level mean field approach can 
only describe the transfer of condensate atoms into the reso- 
nance level <p le!i (r). The subsequent decay of this metastable 
state has been modelled in Ref. [25] in terms of a constant de- 
cay width. Corrections to this approach, which account for the 
energy dependence of this decay width have been provided in 
Refs. 011. 



B. Experimental conditions 

We have solved Eq. (I16> and Eqs. d28t and J29I using a 
cylindrical trap potential with the radial (axial) trap frequen- 
cies of v r = 150 Hz (v z = 1500 Hz) for the sodium experi- 
ments |21], and of v r = 17.5 Hz (v, = 6.8 Hz) in the case 
of 85 Rb 1 2]. As the initial condensate mean field we have 
used the ground state of the Gross-Pitaevskii equation with 
N = 900000 23 Na atoms and N = 1594 85 Rb atoms jH, 
respectively. The initial scattering length was chosen to be the 
background scattering length a\, g for the 23 Na atoms (cf. Ta- 
ble and a = 228flBohi in the case of 85 Rb, corresponding 
to the experimental evaporation field strength of about 162.2 
G. Under these initial conditions the condensate is weakly 
interacting ( 23 Na: « p kfl 3 ~ 6 x 10~ 5 with a peak density of 
n pk a: 1.5 x 10 15 cirT 3 ; 85 Rb: n pk a 3 * 10~ 6 with a peak den- 
sity of « pk « 2.2 x 10 12 cnr 3 ). 

Starting from 162.2 G, in the case of 85 Rb, we have studied 
linear ramps of the magnetic field strength across the 155 G 
Feshbach resonance down to 132 G [2]. In the course of these 
studies we have varied the inverse ramp speeds within the ex- 
perimental range |2] between 1.2/is/G and 200 fis/G. The 
sodium ramps where extended over much larger timescales in 
the experiments 12 ll 15011 . We have therefore chosen the ramp 
intervals in such a way that for each ramp speed the detuning 
\E leti (B)\ of the closed channel resonance state was sufficiently 
large at the initial and final times that increasing the ramp in- 
tervals did not significantly change the total loss of atoms from 
the condensate. 



C. Predictions of the dilferent many-body approaches 

Figure@]shows the theoretical condensate loss curves (N - 
J dx\W(x, t)\ 2 )/N in dependence on the inverse ramp speed for 
the 23 Na Feshbach resonances l5lll . together with the exper- 
imental data from Ref. j^J . In the course of these studies, 
we have solved Eq. dl6l with a variety of different meth- 
ods. These methods involve complete solutions of Eq. J 16b 
including the potential of a spherical trap with the mean fre- 
quency Vho = -^/v 2 v z = 697 Hz, as well as the local density 
approximation, which consists in solving Eq. d!6i in the ho- 
mogeneous gas limit with densities weighted according to the 
initial density distribution associated with the realistic cylin- 
drical atom trap. We have performed these calculations in- 
cluding the two-body time evolution associated with the full 
two channel Hamiltonian © and with the single channel ap- 
proach described in Refs. 1231 12811 . The single channel ap- 
proach includes the time evolution of two atoms in the open 
channel with the background scattering potential adjusted in 
such a way that it recovers, at each magnetic field strength, 
the scattering length in Eq. (I12> as well as the Gribakin and 
Flambaum estimate of the binding energy Et in Eq. J14I . We 
note that the single channel approach does not explicitly in- 
clude the closed channel resonance state. For comparison, we 
have also solved Eqs. J28l > and J29l > of the two level mean 
field approach using a spherical trap potential with the mean 
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♦ spherical trap, 2 channel 

A spherical trap, 2 level mean field 

— cylindrical trap (local density), 2 channel 

— cylindrical trap (local density), 1 channel 

• experimental data 
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Inverse Ramp Speed [[is/G] 



FIG. 4: Loss of atoms from a 23 Na condensate at the end of a lin- 
ear magnetic field ramp across the 853 G and 907 G Feshbach res- 
onances, respectively, as a function of the inverse ramp speed. The 
filled circles are the experimental data from Ref. l2lll . The diamonds 
are determined from Eq. i 1 6i using a spherical trap potential with 
the mean experimental trap frequency. The solid curve represents 
the calculation in the local density approximation, solving Eq. I16i 
for uniform gases, with densities weighted according to the initial 
distribution in a cylindrical trap. The dashed curve shows the cor- 
responding result from a single channel calculation as described in 
Refs. I23T.I28I1 . The triangles correspond to the predictions of the two 
level mean field approach of Eqs. <28> and <29t with a spherical trap 
potential. 



frequency of the cylindrical trap in the experiment 1 2 1 ] . 

The comparison between the theoretical predictions shows 
that although the description of the underlying microscopic 
collision physics differs decisively among the approaches, the 
predicted total loss of condensate atoms is always virtually the 
same. The remnant condensate fractions are therefore largely 
insensitive not only to phenomena related to the inhomogene- 
ity of realistic Bose-Einstein condensates, like collective ex- 
citations, but also to the two channel nature of the two-body 
energy states at magnetic field strengths in the vicinity of the 
Feshbach resonances. In particular, the applicability of the 
single channel approach of Refs. ll23l l28il appears surprising 




50 100 150 

Inverse Ramp Speed [iis/G] 



FIG. 5: Loss of atoms from a 85 Rb condensate at the end of a linear 
ramp of the magnetic field strength across the 155 G Feshbach reso- 
nance as a function of the inverse ramp speed. The circles with error 
bars are the experimental data from Ref. |2]. The squares are deter- 
mined from Eq. <16> using a cylindrical trap potential with the exper- 
imental trap frequencies. The diamonds are obtained with a spheri- 
cal trap potential that accounts only for the mean trap frequency and 
the stars represent a local density approximation for the cylindrical 
trap. The triangles correspond to the two level mean field approach 
of Eqs. 1281 and I29t in a spherical trap. The inset shows the dy- 
namics of the loss of condensate atoms in the cylindrical trap for the 
20 yus/G ramp. 



because the underlying single channel Gribakin and Flam- 
baum estimate of the near resonant binding energy 
in Eq. (I14> is rather poor in the case of the 23 Na Feshbach 
resonances (cf. Fig.|5Jl. The comparison between theory and 
experiment suggests good agreement for the 853 G Feshbach 
resonance of 23 Na but there are significant deviations for the 
broader 907 G resonance. Similar observations were reported 
also in Ref. 12511 . 

Figure|5]shows the analogous theoretical predictions of the 
different approaches with respect to the loss of condensate 
atoms (TV - j dx |T(x, t)\ 2 )/N in a gas of 85 Rb, together with 
the associated experimental data from Ref. Q], in dependence 
on the inverse ramp speed. Also in this case, all many-body 
approaches, including the two level mean field model, yield 
virtually the same results. This is particularly remarkable be- 
cause the physical picture of a single linear curve crossing 
[46], underlying the two level mean field approach, is not ap- 

The 



plicable to the experiments [2] (cf. Subsection IIII Dl 
predicted monotonic curves in Fig.|5]follow the experimental 
trend but they consistently underestimate the observed losses 
in Ref. H. 



VI. UNIVERSALITY OF MOLECULAR DISSOCIATION 
SPECTRA 

Upward ramps of Feshbach resonance levels also play an 
important role in the detection of ultracold molecules pro- 
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FIG. 6: Dissociation spectra of the highest excited vibrational bound 
state of 85 Rb2 as a function of the relative energy of the fragments. 
The speeds of the linear ramps across the 155 G Feshbach resonance 
were varied between 50 G/ms (solid curve) and 1000 G/ms (double 
dashed dotted curve). The ramps started from the experimental evap- 
oration magnetic field strength of 6, = 162.2 G |2] down to a final 
magnetic field strength Bf < Bq = 154.9 G, which was chosen suf- 
ficiently small that at each ramp speed the spectrum was insensitive 
to variations of Bf. The exact calculations based on Eq. 1301 were 
performed with the two channel Hamiltonian J4j. We note that the 
relative kinetic energies of the molecular fragments are given on a 
logarithmic scale. 



duced in atomic Bose-Einstein condensates, as they are fre- 
quently used for the purpose of molecular dissociation upon 
passage across the resonance |y, 0, LL5|, 13811 . The dissocia- 
tion mechanism can be readily understood from the binding 
energies in Fig^for the example of 85 Rb; tuning the mag- 
netic field strength across the resonance position from the high 
field side adiabatically transfers the bound state </>b into the 
zero energy state of two asymptotically free atoms in the open 
channel. Starting from the molecular bound state </>b(B,) at the 
magnetic field strength Z?,, the probability of detecting a free 
pair of atoms with a relative energy between E and E + dE is 
given by Q: 

n(E)dE = p 2 dp j dCl p {(^(BfWiBitfMMBdf ■ (30) 

Here <p${Bf) is the continuum energy state of the two-body 
two channel Hamiltonian @ |23], at the final magnetic field 
strength Bf of the ramp, which is associated with the final 
relative momentum p of the atom pair, and <iQ p denotes the 
angular component of dp. We note that the relative kinetic 
energy of the atoms is related to their relative momentum by 
E = p 2 /m. 

Figure [6] shows the dissociation spectra n{E) as ob- 
tained from Eq. J30l > for the dissociation of 85 Rb2 Feshbach 
molecules in the highest excited multichannel vibrational state 
0b upon passage across the 155 G Feshbach resonance. In 
these calculations the initial magnetic field strength was cho- 
sen as the experimental evaporation magnetic field strength of 



800 



600 



J400 



200- 




400 600 
Ramp speed [G/ms] 



1000 



FIG. 7: Mean single particle kinetic energies of the molecular frag- 
ments associated with the exact dissociation spectra obtained from 
Eq. <30> in Fig. |S| (circles) in dependence on the ramp speed. The 
solid curve has been obtained from the asymptotic prediction in 
Eq. 1321 for comparison. The single particle kinetic energy is half 
the energy of the relative motion of the fragments in Fig. [6] 



Bj = 162.2 G reported in Ref. [2]. We have varied the speeds 
of the linear ramps from 50 G/ms to 1000 G/ms, which il- 
lustrates that the dissociation spectra become broader in the 
energies with increasing ramp speeds. 
The mean kinetic energies 



(£kin) 



1 r 

2 Jo 



dE En(E) 



(31) 



of the molecular fragments after the dissociation characterise 
the speed of expansion of the gas of molecular fragments be- 
fore the detection in related experiments |y, ll5j, 13811 . We 
note that the kinetic energy of a single atom is E^m = p 2 /2m, 
which is half the energy of the relative motion of a pair. A 
Fermi Golden Rule argument 1 15] shows that (Zikin) can be 
expressed in terms of physical parameters of a Feshbach res- 
onance, of the ramp speed \dB/dt\ and of Euler's F function 
MM: 



\£kin) - 



1 



3 H 2 h\dB/dt\ 

4 "V m |flb g l|AB| 



2/3 



r(2/3). 



(32) 



The derivation of Eq. i'32l given in Ref. fl5Tl relies upon 
the linear curve crossing picture described in Refs. I22LI46I1 . 
which identifies the stable multichannel molecular bound state 
</>b with the closed channel resonance state res (r) below the 
dissociation threshold of the background scattering potential. 
Above the threshold les ('") acquires a decay width dependent 
on the energy E KS (B), which is associated with <f> res (r). This 
Fermi Golden Rule decay width leads to the dissociation of 
the molecules. 

We expect Eq. fl32i to be accurate provided that the ini- 
tial and final magnetic field strengths of the linear ramp 
are sufficiently far from the resonance position. This has 
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been demonstrated in the experiments of Ref. fl5ll . in which 
23 Na2 molecules in the highly excited multichannel vibra- 
tional bound state 0b were dissociated upon passage of the 
magnetic field strength across the 907 G Feshbach resonance. 
The derivation in Ref. 11511 is obviously inapplicable to 85 Rb 
ramps in the experimental range of magnetic field strengths 
from 162.2 down to 132 G of Ref. Ql because the reso- 
nance level stays always above the threshold and the molecu- 
lar bound state 0b has only a small admixture of res ('') (< 30 
% in Fig[3Jl. Figure shows, however, that the analytic for- 
mula fl32i is remarkably accurate. Small deviations between 
Eq. fl32i and the exact mean energies obtained from Eqs. J30I 
and OH only occur with increasing ramp speeds. The molec- 
ular dissociation spectra thus exhibit the same insensitivity to 
the exact microscopic binary collision physics as the remnant 
condensate fractions in Section IV1 

VII. CONCLUSIONS 

We have shown in Sections [V] and I VII that remnant con- 
densate fractions or molecular dissociation spectra associated 
with linear upward ramps of a Feshbach resonance level do 
not reveal the details of the underlying microscopic binary 
collision physics. The formula fl32i for the mean single par- 
ticle kinetic energy of molecular fragments suggests that the 
Landau Zener parameter I22H23I1 

N-m \a bg \\AB\ 

olz = 47r (33) 

<V m \dB/dt\ 

determines the universal observables in Sections [V] and IV1I 
when |flbgL I AZ?| and \dB/dt\ are the parameters that can be 
varied independently. Here N is the number of atoms and *V 
is the volume of the gas. For N - 2 atoms the Landau Zener 
coefficient exactly determines the final population in the zero 
energy level by 

po = exp(-27r<5 L z) (34) 

in a two level linear curve crossing model of a Feshbach res- 
onance i23ll46Tl . provided that the linear ramp starts and ends 
asymptotically far from the crossing point B res (see Fig. [Q. 
The analysis of the two level mean field equations d28l and 
d29l of Ref. ll23ll reveals that the final populations of the rem- 
nant condensate and of the closed channel resonance state are 
indeed determined by 6lz, i e- the final populations are con- 
stant when the parameters |abgL I AZ?| and \dB/dt\ are varied in 
such a way that Slz is kept constant. In particular, the depen- 
dence of Eq. fl33i on the modulus of the ramp speed shows that 
the final populations are insensitive with respect to the ramp 
direction. 

The analysis in Ref. ll52il shows that Eq. J34I determines 
the exact asymptotic population po in the open channel ground 
state even when the complete quasi continuum of modes in the 
open channel is taken into account in a linear curve crossing 
model of a Feshbach resonance. This explains why including 
the full two channel time evolution operator in Eq. dl6l leads 
to the same results as the two level mean field approach for 



the 23 Na Feshbach resonances in Section IV1 given that phe- 
nomena related to the non Markovian nature of Eq. dl 6t are 
negligible. Our studies of the 155 G Feshbach resonance of 
85 Rb in Figs.|5]and0suggest that also the asymptotic criterion 
of applicability of the two-body Landau Zener approach can 
be violated to a large extent without any significant changes 
in the results. We note, however, that the two level pic- 
ture is not suited to describe the interferometric experiments 
of Refs. ||3l 0], which employ a sequence of magnetic field 
pulses. In fact, the occupation of binary continuum modes, 
observed in these experiments, is ruled out in any two level ap- 
proach from the outset and, in addition, the nature of the high- 
est excited vibrational diatomic bound state is not recovered. 
This explains why the two level mean field approach fails to 
predict the experimentally observed atom-molecule Ramsey 
fringes 13911 - 

Although the Landau Zener coefficient in Eq. d33i can be 
derived from a linear curve crossing model of a Feshbach res- 
onance, it is independent of the slope dE res /dB of the res- 
onance and contains just the parameters flb e and AB of the 
Feshbach resonance, which determine the scattering length in 
Eq. H2\ . In fact, the inset of Fig. |5] shows that the main loss 
of condensate atoms occurs in the close vicinity of the Fes- 
hbach resonance, where all low energy scattering properties 
become universal |23, 43], i.e. the associated physical observ- 
ables just depend on the scattering length a. The universal 
parameters «bg and AB as well as the inverse ramp speed that 
determine Eq. (1331 are also accounted for in the single channel 
description of the first order microscopic quantum dynamics 
approach of Ref. 12811 . We have, indeed, confirmed that the 
single channel description also recovers the theoretical loss 
curves in Fig. [5] as it does in the case of the 23 Na Feshbach 
resonances in Fig. @] Consequently, our studies clearly re- 
veal that the losses of condensate atoms in linear ramps of the 
magnetic field strength across Feshbach resonances depend, 
to a large extent, just on the Landau Zener parameter. The 
two or multichannel nature of Feshbach resonances is virtu- 
ally irrelevant to the loss curves in Figs. 0] and |5] 

In summary, we have shown that physical observables, like 
the loss of condensate atoms or molecular dissociation spec- 
tra, associated with linear upward ramps of a Feshbach res- 
onance level, are largely independent of the microscopic bi- 
nary collision dynamics and depend just on the parameter 
\a\ )% \\AB\l\dB I dt\ and on the local densities of the atoms in 
the gas. Consequently, related experimental studies are in- 
conclusive with respect to the details of the interatomic inter- 
actions. The interferometric studies in Refs. |3j|41 and their 
subsequent theoretical analysis in Refs. [26, 28, 321 13911 have 
demonstrated, however, that the universality of the physical 
quantities described in this paper is characteristic for linear 
ramps of the magnetic field strength. 
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